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Abstract. The Coulomb gauge Hamiltonian model is used to calculate masses for selected J PC states 
consisting of exotic combinations of quarks and gluons: ggg glueballs (oddballs), qqg hybrid mesons and 
qqqq tetraquark systems. An odderon Regge trajectory is computed for the J glueballs with intercept 
much smaller than the pomeron, explaining its nonobservation. The lowest 1~ + hybrid meson mass is found 
to be just above 2.2 GeV while the lightest tetraquark state mass with these exotic quantum numbers is 
predicted around 1.4 GeV consistent with the observed 7r(1400). 

PACS. 12.38.Lg Other nonperturbative calculations - 12.39.Ki Relativistic quark model - 12.39.Mk 
Glueball and nonstandard multi-quark/gluon states - 12.40.Yx Hadron mass models and calculations 



1 Introduction 

Establishing the existence of exotic hadrons (non qq or 
qqq structure) is of paramount importance and remains 
one of the key unsolved problems in hadronic physics. In 
particular, it is expected from general QCD principles that 
nonconventional color singlet states of gluons and quarks 
should exist such as glueballs {gg and ggg), hybrid mesons 
{qqg) and multiquark states {qqqq and qqqqq). The present 
work addresses the structure of these hadrons and provides 
new information to assist experimental searches. 



2 Coulomb Gauge Hamiltonian Model 

In a series of publications ^ni21I3l^3l5lEll3IHI a realistic model 
for hadron structure has been developed and applied to 
the quark and gluon sectors. This field theoretical, rel- 
ativistic many-body approach utilizes an effective QCD 
Hamiltonian, H e s, formulated in the Coulomb gauge. It 
properly incorporates chiral symmetry using standard bare 
current quark masses but dynamically generates a con- 
stituent mass and spontaneous chiral symmetry breaking 
2 . Through approximate many-body diagonalizations it 
successfully describes the meson spectrum ^JE] and is con- 
sistent [7] with lattice glueball predictions. It also yields 
a good description of the vacuum properties (quark and 
gluon condensates) within a minimal two parameter the- 
ory. 

The effective Hamiltonian, an approximation to the 
exact Coulomb gauge QCD Hamiltonian, is 



H q + H g + H qg + He 



H q = J dxW^ (x)[- ia ■ V + (3m]ip(x) 



(1) 
(2) 



H g = \ Jdx [IF (x) • n a (x) + B Q (x) • B a (x)] (3) 
H qg =gJdxJ a { X )-A a {x) (4) 



Hc = -^\ dxdyp a (x)V{\x~y\)p a {y) 



(5) 



with g the QCD coupling, W the quark held, m the current 
quark mass, A a the gluon holds satisfying the transverse 
gauge condition, V • A a = 0, a = 1, 2, ...8, 11° the conju- 
gate fields and B a the non-abelian magnetic fields 



B a = V x A a + ig/ Qbc A b x A c . 

The color densities, /? a (x), and currents, J a , are 
p Q (x) = «^(x)T°<F(x) + / Qbc A b (x) • n c (x) 



J a = tf rt (x)aT a tf'(x) 



(G) 

(7) 
(8) 



with T a — 4p and f abc the SU3 color matrices and struc- 
ture constants, respectively. Confinement is described by 
a Cornell type potential, 

V(r=\x-y\) = V c (r) + V L (r) (9) 

V c (r) = (10) 
r 

V L {r) = or, (11) 

with previously determined string tension, a = 0.135 GeV 2 , 

2 

and a s = j- = 0.4. Below we denote the Fourier trans- 
form of V by V. 

The bare fields have the Fock operator expansions (quark 
spinors u, v, helicity, A = ±1, and color vectors £0=1.2,3) 



dk 

"2^)3 



!ffc(k)e lk - x e c 



(12) 



2 
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A Q (x) 



WA (k)& AC (k)+ UA (-k)4 c (-k) 
(^^f [aa(k) + aat( - k)]eIk - : 



na(x) - z /(^/l [aQ(k) - aat( - k)]eI 



(13) 
(14) 

(15) 



with quark, anti-quark and gluon Fock operators b\c(k), 
d\c{— k) and a^(k) (/i = 0, ±1), respectively. The Coulomb 
gauge condition, k-a a (k) = (— l^fc^a™ (k) = 0, produces 
transverse commutation relations, 



3 Applications 

Predictions for the low- lying spectra of glueballs, hybrid 
mesons and tetraquark systems are now presented and 
discussed. Since these hadrons consist of 3 or more con- 
stituents, the masses for selected J PC states are computed 
variationally 



Mjpc 



(V JPC \H cS \y JPC ) 

(tyJPC\tyJPC} 



(23) 



a«(k),<#(k')] = (2^) 3 <U<5 3 (k- k')IV(k) , (16) 



with 



D fifi ,(k) = S u 



V ' 7,2 



(17) 



The ground state (model vacuum) is generated us- 
ing the Bardeen-Cooper-Schriffer (BCS) method, entail- 
ing rotated field operators (Bogoliubov-Valatin transfor- 
mation), 

S AC (k) = cosy6 A c(k) - Asinydj^-k) 
D A c(-k) = cosyd AC (-k) + Asin yb^ c (k) 

Pi, Q a t 



The variational approximation has been comprehensively 
tested in two body systems by comparison with exact di- 
agonalization and found to be accurate to a few percent. 



3.1 Glueballs 

Previous work jS] has investigated gg glueballs which only 
have C = 1. For C = — 1 glueballs (oddballs), Fock states 
with at least 3 gluons are necessary and the variational 
wavefunction is (q;=i,2.3 are the cm gluon momenta) 



JPC\ 



J dqidq 2 dq 3 <5(qi + q2 + q 3 ) 



(24) 



a°(k) = coshe fc a"(k) +sinh0 fe a^(- k ) , (18) ^^(ca^.^C^Q-t^J^toJa^Cq,)!^ 



producing the dressed, quasi-particle operators a a , B\c 
and D\c, respectively. The quasi-particle (BCS) vacuum, 
determined by B\ C \0) = D X cW) = ofyW) = 0, is built 
on the bare parton one, 6 A c|0) = d\c\0) = a^\0) = 0, 

\Q quark ) = e -/l^ Atan * b ^(k)4 c (-k) |0) (ig) 
tanhe h D^(k)o2t(k)o"I(-k), 



= p_ J (2 



I ^ gluon ) G 

The composite BCS vacuum, |J?) = \f2 quar k 



10). (20) 



I ^ gluon ) j 

contains quark and gluon condensates (correlated qq and 
gg Cooper pairs). A variational minimization of the vac- 
uum expectation value of the Hamiltonian, 5{f2\H c s\fi) — 
0, yields the constituent quark and gluon gap equations 

2 f_dq_ 

3 J (2nY 



ks k - mc k 



-(skCqX - s q c k )V{\k- q|) (21) 



fc 2 
3 



dq 

(2^3 
dq 1 



n|k-q|)[l 



( 22 ) 

(27T) J LU q 

with Sfc = sin4>k, = cosipk and a; = k-q. Here 0j. = <^>(fc) 
is the quark gap angle related to the BCS angle by, 
tan(4>k — 0k) = m/k, and ujk = ke~ 2Bk is the effective 
gluon self energy. The last term in Eq. I122H is due to the 
non-abelian component of the gluon kinetic energy. The 
quark gap equation is UV finite for the linear potential 
since Vt(|p|) = Sna/p 4 , but not for the Coulomb po- 
tential Vc(|p|) = — 47ra s /p 2 . The gluon gap equation has 
both logarithmical and quadratical UV divergences and 
an integration cutoff, A — 4 GeV, determined in previous 
studies is used in both equations. 



with color tensor C abc either totally antisymmetric f abc 
(for C = 1) or symmetric d ahc (for C = —1). Boson 
statistics thus requires the C = — 1 oddballs to have a 
symmetric space-spin wavefunction. Using eq. and a 
two-parameter variational radial wavefunction, the J 
oddball states have been calculated. Only the Abelian 
component of the magnetic fields ae retained and the hy- 
perfine interaction, H qg , is suppressed. There are three 
terms contributing to the mass expectation value which 
are depicted in Fig. ^ which correspond to the gluon self- 
energy (top), gluon-gluon scattering (middle) and annihi- 
lation (bottom). 

The nine-dimensional variational calculation was per- 
formed using the Monte Carlo method with the adap- 
tive sampling algorithm VEGAS !) and numerical conver- 
gence required between 10 5 and 10 6 samples. The oddball 
mass predictions are compared in Tabled to available 
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Fig. 1. Glueball diagrams for (^ g p g c \H cS \^ g J g p g c ). 
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Table 1. Glueball quantum numbers and masses in MeV. 
Error in H c g (from Monte Carlo only) is 100 MeV or less, the 
quoted lattice errors are typically 200-300 MeV. 



Model 


1" 3"- 5— 7— 


Coulomb gauge H c g 
lattice [101 
lattice [TT] 
Wilson-loop [H 


3950 4150 5050 5900 
3850 4130 
3100 4150 
3490 4030 



lattice gauge results |10llll| and a Wilson-loop inspired 
model A study of the glueball mass sensitivity to 
both statistical and variational uncertainties yielded error 
bars at the few per cent level. 

Figure displays predicted oddball Regge trajectories 
from the alternative approaches. Lattice results are de- 
picted by open circles JOj and diamonds ^J. Constituent 
gluon predictions are represented by boxes, solid trian- 
gles and solid circles and correspond to a Wilson-loop in- 
spired potential model |12|. a simpler harmonic oscillator 
calculation labeled Hm and the H e s approach, respec- 
tively. The odderon trajectories for the latter two models 
are represented by the solid lines, Qq = 0.18i + 0.25 and 
a< = 0-23it — 0.88, which provide an overall theoretical 
uncertainty. The much steeper dashed line is the u> trajec- 
tory. 



3.2 Hybrid mesons 

We denote the momenta of the dressed quark, anti-quark 
and gluon by q, q and g, respectively, and work in the 
hybrid cm system. The color structure for a qqg hybrid is 
given by SU C {3) algebra, (3 ® 3) ® 8 = (8 <g> 8) © (8 <g> 1) = 
27©10©10©8©8ffi8ffil. Hence for an overall color singlet 
the quarks must be in an octet state like the gluon which 
leads to a repulsive qq interaction, confirmed by lattice at 
short range, that raises the mass of the hybrid meson. The 
hybrid wavefunction has the general form 

\*ifg ) = /dqdqdg *(q + q + g) *& C (q, q, g) 

T- s B{ c ( q )Dl s m)a^( S )\n) . (25) 

We have extended our previous hybrid study |2] by in- 
cluding the H qg Hamiltonian term containing the J a • A a 
operators. Following jS], an effective quark hyperfine inter- 
action with a J° • J° form is obtained using perturbation 
theory to second order in g and integrating over the glu- 
onic degrees of freedom. This contribution to the hybrid 
mass is represented by the qq gluon exchange Feynman 
diagrams in Fig. |31 (first two in the bottom row). The non- 
abelian magnetic field terms are also included and entail 
triple-gluon vertices (last two diagrams in Fig. [3} . The re- 
maining diagrams represent the self-energy, scattering and 
quark annihilation mass contributions. 



Oddball Chew-Frautschi plot 




Morningstar and Peardon 
Meyer and Teper 
Kaidalov and Simonov 
(O trajectory 
H M (this work) 
Coulomb gauge H (this work) 



M 2 (GeV 2 ) 

Fig. 2. Odderon trajectories from constituent gluon models 
and lattice compared to the ui meson Regge trajectory. 



Three key results follow. First, the predicted odderon 
has slope similar to the pomeron but intercept clearly 
lower than the to value. Second, the odderon starts with 
the 3~~ state and not the 1~~ which is on a daughter 
trajectory. Note that there are no lattice 5 glueball pre- 
dictions which are necessary to confirm this point and we 
strongly recommend that future studies calculate higher 
states. Third, all approaches agree that the 3~~ mass 
is near 4 GeV. 
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Fig. 3. Hybrid meson diagrams for (^^ c \H e s \&qq g 



The hyperfine interaction from the H qg term is 

v t = \J jdxdyJ? W^(x,y)J°(y), 
with kernel reflecting the transverse gauge 



U(\x-y\ 



(26) 



(27) 



The potential U is a modified Yukawa with dynamical 
mass, m g = 600 MeV, for the exchanged gluon as ex- 
plained in 0. Its momentum space representation is 



U(p) 



8.04 
7? 



2 (-^ T +0.82) 



V ln°- 8 (-£l + 1.41) P>m 9 _ (28) 



24.50 
p 2 



P < m g 
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The quark hyperfine interaction also generates additional 
terms in the quark gap equation 0. 

The 12 dimensional integrals were calculated using the 
Monte Carlo method and repetitively evaluated with an 
increasing number of points until a weight-averaged result 
converged, typically involving about 50 million samples. 
The hybrid mass error introduced by this procedure is 
about ± 50 MeV. For each J PC hybrid state we optimized 
the two variational parameters. 

Using standard current quark masses, m u — rrid = 5 
MeV, m s = 80 MeV and m c = 1000 MeV, the pre- 
dicted low-lying mass spectra for light and heavy hybrid 
mesons are presented in Figs. 01 and El respectively. Note 
that quark annihilation interactions increase the hybrid 
mass and this introduces isospin splitting since it only 
contributes in the I q q channel. More importantly, all hy- 
brid masses, especially the lightest exotic 1 h state, are 
clearly above 2 GeV. This is consistent with lattice jl.'il 
raiODl and Flux Tube model I5T1I2^II2^I re- 



sults summarized in Table El These composite predictions 
strongly suggest that observed 1 h 7r(1600), and more 
clearly 7r(1400), are not hybrid meson states. 



3.0 



2.5 



2.0 



^exotica 


^exotica 


3" + 






3~ + 


- — r + 




— 2 ++ — r + 


— r + 
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uOg (1=1) 



Fig. 4. Low lying uug spectra. 



M(GeV) 
3.0 

2.5 

2.0 



exotica 



4.5 



-l; 4.o 
-i 

-o" 



3.5 



exotica 



-o~x 
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Fig. 5. Low lying ssg and ccg spectra. 

The charmed ccg hybrid spectrum has a slightly dif- 
ferent level order compared to the strange ssg and uug 
spectra due to the hyperhne interaction. The predicted 
strange and charmed exotic 1 ^ states are also in reason- 
able agreement with both lattice and Flux Tube results. 



Table 2. Published predicted exotic 1 + masses, in GeV, for 
light, strange and charmed hybrid mesons. 



Model 


u/d hybrid 


s hybrid 


c hybrid 


Lattice QCD [13-20] 
Flux Tube 21 22|23| 


1.7 - 2.1 

1.8 - 2.1 


1.9 
2.1 - 2.3 


4.2 - 4.4 
4.1 - 4.5 



3.3 Tetraquark systems 

This is the first four-body application using this approach. 
The SU C (3) color algebra for four quarks produces 81 color 
states, 3®3®3®3 = 27©10ffil0ffi8ffi8ffi8ffi8©lffil, of 
which two are color singlets that can be obtained in four 
different ways, depending on the intermediate color cou- 
pling: singlet scheme (molecule), octet scheme, and two 
diquark schemes involving the triplet and the sextet rep- 
resentations (see Fig. El . 



q 



q 


q 


2 


q 



8 \ ® 8 = 27®10©10ffi8®8®l 
] ® \ l /= l (molecule) 



/6\® /6V= 27 I 8 • 1 

I $ | i ffi 

I 3 I ® I 3yt= 8*1 



Fig. 6. Color singlets via four different representations. 

Working in the cm system and denoting the momenta 
of the quarks by qi and q3 , and those of the anti-quarks 
by q2 and q4, the tetraquark wavefunction is 



K PC ) 



dqidq 2 dq 3 dq 4 (5(qi + q 2 + q 3 + q4) 



IPC 

^i A2A3 Az 



(q 2 )S 



A 3 C 2 



(qi,q2,qs,q4)i?c3cl 

{^s)D{ iCi {cu,)\fi quark ) 



(29) 



where the color elements Rc 3 cl depend on the specific 
color scheme chosen. 

Contributions to the Hamiltonian expectation value 
are summarized in Fig. [7| and correspond to 4 self-energy, 
6 scattering, 4 annihilation and 70 exchange terms each 
of which can be reduced to 12 dimensional integrals that 
are evaluated in momentum space. Because of the com- 
putationally intensive nature of this analysis, the hyper- 
fine interaction was not included. Performing large-scale 
Monte Carlo calculations (typically 50 million samples), 
has conclusively determined that the molecular represen- 
tation (i.e. meson-meson) produces the lightest mass state 
for a given J PC . This is due to suppression of certain inter- 
actions (cancellation of color factors) in the singlet-singlet 
molecular representation and also the presence of repul- 
sive forces in the other, more exotic, color schemes. 

Using m u = md = 5 MeV, the predicted tetraquark 
ground state is the non-exotic vector 1 ++ state in the 
molecular representation with mass around 1.2 GeV. Fig- 
ures El and El depict the predicted tetraquark spectra for 
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Fig. 7. Tetraquark diagrams for (W&f ° \H eS ° ) 
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Fig. 8. Tetraquark 1=1 singlet scheme (molecule) spectra. 



states having conventional and exotic quantum numbers 
in both singlet and octet color representations. Due to 
quark annihilation interactions (qq — > <? — > qq) in the 
Iqq = channel there are isospin splitting contributions, 
up to several hundred MeV, in the octet but not singlet 
scheme as illustrated in the figures. The annihilation inter- 
action terms are repulsion, yielding octet states with 1 = 2 
lower than the 1 = 1 which are lower than the 1 = 0. The 
molecular-type states are all isospin degenerate and the 
lightest exotic molecule is the with mass 1.35 GeV. 
Because of the isospin degeneracy, there will be several 



molecular-type tetraquark states with the same J in 
the 1 to 2 GeV region. Further, these states can be ob- 
served in different electric charge channels (different I z ) 
at about the same energy, which is a useful experimen- 
tal signature. There are 3 orbital angular momentum and 
the lightest 1 h exotics have only 1 p-wave. The lightest 
1 h is predicted near 1.4 GeV which is close to the ob- 
served 7r(1400), suggesting this state has a meson-meson 
type structure distinct from states having quarks in more 
exotic color representations (quark atoms). Related, the 
computed tetraquark mass for 1 ^ states with these more 
exotic color configurations is above 2 GeV. This is consis- 
tent with model predictions [S] for exotic hybrid meson 
(qqg) 1~ + states also lying above 2 GeV due to repul- 
sive color octet quark interactions. Finally, for the same 
J PC states, including the 1 the computed masses (not 
shown) in both the triplet and the sextet diquark color 
representations are all heavier than in the singlet repre- 
sentation and comparable to the octet scheme results. 



M(GeV) 
2.5 — 



2.0 



1.5 



-1 

"0 



lO ,2 



; Octet color scheme I 



1=1 (1x1) 1=1 (1x0) 

Fig. 9. Tetraquark 1=1 octet scheme spectra. 



4 Conclusions 



Summarizing, using the Coulomb gauge Hamiltonian model 
we have performed large-scale Monte Carlo calculations 
for ggg glueballs, qqg hybrid mesons and qqqq tetraquark 
systems. Our oddball spectrum, along with lattice and 
other glueball model predictions, clearly documents an 
odderon trajectory with slope similar to the pomeron but 
much lower intercept. This explains why the odderon has 
not been seen in reactions with pomeron exchange. If the 
odderon intercept is comparable to the u> value it may be 
possible to observe it in reactions where the pomeron is 
absent such as pseudoscalar [21] or tensor meson [2H| elec- 
tromagnetic production. However if our prediction that 
the intercept is below 0.5 is correct, it is unlikely the odd- 
eron will be seen. It is important therefore that lattice 
calculations for the 5 be performed to confirm this. 

For the hybrid meson investigation, the Coulomb gauge 
Hamiltonian model predicts that all uug and udg states 
have mass above 2 GeV. In particular, lattice and con- 
stituent model predictions for light, strange and charmed 
exotic 1 hybrid mesons are collectively in agreement. 
This strongly suggest that the 7Ti(1400) and 7Ti(1600) are 
not hybrid mesons but have an alternative structure as 
argued below. 

Our tetraquark results clearly show that [(qq)i®(qq)i]i 
molecular states are lighter than the more exotic color 
octet [(qq) g ® (qq) s]i atomic-like states. Most significantly, 
for the 1 h , / = 1 channel our predicted lightest color 
octet state is above 2 GeV, while the lightest singlet scheme 
state is around 1400 MeV, close to the 7r(1400). Combining 
these results with several model 1 h hybrid meson predic- 
tions indicates that the observed ir states are not hybrids 
or exotic color configurations of 4 quarks but rather more 
conventional meson-meson type molecules. 

Future work will address charmed (ccuu) and strange 
tetraquark (cess) systems to study the recently reported 
X(3872) and V(4260) states. Dynamical mixing of glue- 
ball, hybrid and conventional meson and tetraquark states 
will also be investigated along with three-body forces [2E] ■ 

The authors are very grateful to F. J. Llanes-Estrada and P. 
Bicudo for insightful discussions. Work supported by U. S. 
DOE grants DE-FG02-97ER41048 and DE-FG02-03ER41260. 
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